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A plasma transport theory that spans weak to strong coupling is developed from a binary colli- 
sion picture, but where the interaction potential is taken to be an effective potential that includes 
correlation effects and screening self-consistently. This physically motivated approach provides a 
practical model for evaluating transport coefficients across coupling regimes. The theory is shown 
to compare well with classical molecular dynamics simulations of temperature relaxation in electron- 
ion plasmas, as well as simulations and experiments of self-diffusion in one component plasmas. The 
approach is versatile and can be applied to other transport coefficients as well. 
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The microscopic dynamics of Coulomb collisions deter- 
mines macroscopic transport properties of plasmas such 
as difFusivity, resistivity, viscosity, etc. [U 1^. Usually 
plasmas are so hot and dilute that the average particle ki- 
netic energy greatly exceeds the potential energy of inter- 
action. In this weakly coupled regime, Coulomb collisions 
consist of a series of many small angle binary scattering 
events [SHS]- Strongly coupled plasmas are fundamen- 
tally different. In this regime, the interaction potential 
energy exceeds the particle kinetic energies, so scatter- 
ing angles are large and correlation effects are impor- 
tant. Plasmas in several modern experiments, including 
inertial confinement fusion (ICF) [315], antimatter plas- 
mas ultra cold plasmas [TU] and dusty plasmas [TT] . 
exhibit strong coupling effects; as do some naturally oc- 
curring objects including neutron star crusts [12], white 
dwarf stars [TS] [M] and giant planet interiors [T5l - fT7] . 

Understanding how transport properties are modified 
in strongly coupled plasmas is interesting both from a 
fundamental physics standpoint and as a practical mat- 
ter. Accounting for correlation effects remains a chal- 
lenge for theory, even though accurate transport coeffi- 
cients are critical input to the macroscopic (ffuid) equa- 
tions used to model these systems. Transport calcula- 
tions typically rely on computationally expensive particle 
simulations, such as molecular dynamics (MD) [T5l420j . 
Analytic theory is desirable because it can both eluci- 
date the physical processes that inffuence transport at 
strong coupling, and provide an efficient means for es- 
timating the transport coefficients that fluid equations 
require as input j2TU23| . In this Letter, we describe a 
physically motivated method of extending conventional 
transport calculations, which is efficient enough to be 
practically implemented in fluid simulation codes. The 
theory provides coefficients that agree with experimen- 
tal [24| and classical MD simulation data [25l |26] across 
coupling regimes. 

Like weakly coupled theories, our theory is based on 
a binary collision picture, but where particles interact 
via an effective potential that includes average effects 



of the intervening medium; including both correlations 
and screening. This effective potential is used to derive 
a scattering cross section, which is then applied to the 
Boltzmann collision operator and Chapman-Enskog col- 
lision integrals [T] to calculate the various transport co- 
efficients. 

In fact, traditional plasma theories also rely on an ef- 
fective potential. The bare Coulomb potential neglects 
screening of the intervening medium and leads to a di- 
vergent collision operator. To fix this unphysical diver- 
gence. Landau utilized the weak coupling assumption 
to cutoff the impact parameter at the Debye screening 
length Hence, imposing an effective potential. This 
approximation leads to the traditional Coulomb loga- 
rithm. In A, where A ~ r~'^/^ is the plasma parameter. 
It is valid in the limit that this parameter is asymp- 
totically large. Here, the Coulomb coupling parame- 
ter, r = Z'^e'^/ikBTa) where a = {4Trn/3)-^/^ is the 
Wigner-Seitz radius, will be used to quantify coupling 
strength. The Lenard-Balescu equation is an alternative 
plasma kinetic theory derived from the BBGKY hierar- 
chy [^. It has the advantage of accounting for screen- 
ing self-consistently, but it does not account for close 
interactions and also diverges. Again, this divergence is 
fixed through the weak coupling approximation. Using 
the screened Coulomb (Yukawa) potential as an effective 
potential avoids the logarithmically divergent integrals 
and can extend the binary collision approach to larger 
coupling strength [S7HST] . but it does not capture corre- 
lation effects, which onset when L > 1. Can the binary 
collision picture be extended further by using an effective 
interaction potential that accounts for correlation effects 
in addition to screening? In this Letter, we present evi- 
dence that it can. 

Previous theories of transport in strongly coupled plas- 
mas have largely focused on developing new closure 
schemes of the BBGKY hierarchy that include correla- 
tions 32 36 . These typically either derive a new colli- 
sion operator that has a generalized linear dielectric re- 
sponse with local field corrections [35J [33] , or calculate 
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transport properties from higher-order equihbrium cor- 
relation functions [35l |36] . The sahent feature of these 
closures is that they go beyond the mean field approxi- 
mation of conventional plasma theories to include corre- 
lation effects. The mean field approximation is justified 
in the weakly coupled limit because there are many par- 
ticles within the interaction length scale (Debye length) 
of a test particle. However, in a strongly coupled plasma 
the interaction distance is instead characterized by the 
inter-particle spacing. In this regime, the test particle 
self interaction must be neglected, and correlations ac- 
counted for. 

Next, we establish a relationship between the effec- 
tive interaction potential {(f)) and pair correlation func- 
tion (32)- This enables determination of (j) from closures 
that include correlation effects in 92- To illustrate this 
point, consider the second BBGKY equation for 32 for a 
classical one component system 
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and 



v{\ri ~ Tjl) is the bare Coulomb potential. The usual 
kinetic theories can be obtained by neglecting certain 
terms in Ec^. (1). The Landau collision operator is ob- 



tained by neglecting terms ( lb |-( Id I, the Lenard-Balescu 
collision operator by neglecting (lb I and (Ic), and the 
Boltzmann collision operator by neglecting ( Ic ) and ( Id I. 
Each choice defines an approximation for 32 and, in turn, 
a different collision operator. 

For our purposes, the equilibrium limit of these ap- 
proximations is instructive. At equilibrium, 172(1,2) — 
"-^/M(Pi)/M(P2)/i(|ri - r2|) where /m is a Maxwellian, 
h{r) = g(r) — 1 and g{r) is the pair distribution function; 
ng{r) is the average density at a distance |r| from any par- 
ticle. In this limit, the Landau and Lenard-Balescu clo- 
sures correspond to assuming (;'"(r) = 1 — ei;(r)/fcBr and 
g^^{r) = 1 - e4>scir)/kBT where 0sc = gexp(-r/A_D)/r 
is the screened Coulomb potential. These can be ob- 
tained from the weakly coupled limit (ecfi/k-QT <^ 1) of 
the general equilibrium relationship |37| 



where 



V0: 



g{r) = exp{~e(l)/kBT), 
/e-^/fcB^(-V^[/)rfr3...drjv 



(2) 



(3) 



defines an effective interaction potential [38], and U = 
Si j Vi,3 is the total interaction energy. The quantity 
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FIG. 1. (a) Pair distribution function for the OCP deter- 
mined from classical MD (circles), HNC (solid lines) and the 
screened Coulomb potential (dashed lines) for F = 0.1, 1, 10 
and 100. (b) Magnitude of the effective interaction potential 
calculated from Eq. The thick red and blue lines show 
the r/a intervals where <j) <Q for F = 10 and 100. 



represents the mean force acting on particle 1, with 
particles 1 and 2 held at fixed positions (ri and T2), aver- 
aged over the positions of all other particles. The usual 
plasma theories rely on the e4>/kB,T <^ 1 assumption, 
whereas the Boltzmann collision operator does not. Ac- 
tually, the Boltzmann collision operator corresponds to 
using the bare Coulomb potential (0 = v) in Eq. (|2|, 
which neglects both screening and correlations. We will 
use an effective potential that includes these effects. 

Figure [T^ shows a comparison of g{r) obtained 
from different approximations with that extracted from 
classical MD simulations of a one component plasma 
(OCP) [53]. Figure [1^ shows the corresponding effective 
potential from Eq. T2]). The screened Coulomb poten- 
tial is an excellent approximation in the weakly coupled 
regime, F ^ 1, but this breaks down as correlation effects 
onset at F > 1. It fails entirely in the strongly correlated 
cases. To obtain an analytic approximation for g(r) that 
includes correlation effects, but doesn't rely on compu- 
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tationally expensive MD simulations, we use the hyper- 
netted chain (HNC) closure. HNC is a well-established 
approximation in which g(r) is determined from the two 
coupled equations [33] 



5(r) 
and 



exp 



v{r)/kBT + n c{\r-r'\)h{r')dr' (4) 



/i(k) = c(k)[l-Hn/i(k)]. (5) 



For the OCP, v{r)/kBT = Ta/r and ^(k) denotes the 
Fourier transform of h{r). Figure [I] shows that, like 
the screened Coulomb potential, HNC provides an ex- 
cellent approximation for weak coupling. However, un- 
like the screened Coulomb, this agreement extends to the 
strongly coupled regime. Next, we apply this effective 
potential to determine transport properties. 

In the Chapman-Enskog theory, transport coefficients 
arise through the collision integrals |lj 



/ d^^^''+^e~«' / d9ass' sm0{l-cos' 0] 
Jo Jq 

(6) 



in which 9 is the scattering angle, ^ — u/vg 
"^Ts + "^Ts'' '^Ts ~ "^Ts/ras, u — |v — v'l and ffgs' is the 
differential scattering cross section. Here, s and s' denote 
species. Alternatively, these can be written [SU] 



^ 3 rUs Psa' ^gs' 

16 m.,.,' n.,' S.,.,' 



where 



1 

2 7o 



(7) 



(8) 



is a "generalized Coulomb logarithm" associated with the 
(Z, collision integral. Here, S^,' = S^^;^'' is the lowest 
order term, 



is a reference collision frequency. 



POO 
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(10) 



is the l*"^^ momentum-transfer cross section, CTo — 
{nq^q^,) / [mlgivlg,) is a reference cross section, and 
rriss' = Tngmgi / {rrig -\- mg') is the reduced mass. The 
scattering angle is that of a classical binary collision 

e = b drr-^[l-b^/r^-2e(l){r)/{mgg,u^)]-^/^ (11) 



in which Vf, is the distance of closest approach, de- 
termined from the largest root of the denominator in 
Eq. lul). 
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FIG. 2. Self diffusion coefficient for the OCP calculated using 
classical MD with Green-Kubo relations (blue circles), and 
using various effective potentials in the Chapman-Enskog col- 
lision integrals: from MD derived g{r) data (black squares), 
HNC (red diamonds) and screened Coulomb (black dashed 
line). 



Equations ([7l)-(|Tl|) determine the transport coeffi- 
cients. First, we recover the weakly coupled limit. Ap- 
plying the screened Coulomb potential provides gen- 
eralized Coulomb logarithms that avoid the tradi- 
tional divergences pTHHU] . For weak to moderate cou- 
pling, the lowest order Coulomb logarithm is S'^'^^' = 
exp(A~-'^)i?i(A~i^), where Ei is the exponential inte- 
gral [30]. The weak coupling limit of this returns the 
conventional Coulomb logarithm including an order unity 
correction 



^(1,1) 



lnA-7 = ln(0.56A). 



(12) 



This order unity correction extends the conventional In A 
solution to the moderate coupling regime In A > 2. It has 
also been obtained by others using complicated renor- 
malization techniques [391 - I4T] . The result that is usually 
cited is ln(0.765A). Equation (IT2| reduces to ln(0.79A) 



in the limit — Ti and A is defined using the electron 
Debye length, which is within ^ 3% of the coefficient 
from previous calculations. Next, we extend these cal- 
culations into the strong coupling regime using the HNC 
obtained effective potential to compute the self diffusion 
coefficient for a OCP and the temperature relaxation rate 
of an electron-ion plasma. 

Self diffusion in a OCP: The Chapman-Enskog self 
diffusion coefficient to first order is [T] 



[Dg 



16. 



(1,1)' 



(13) 



Accounting for a second order correction resulting from 
deviations from Maxwellian distributions provides 
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l/(1-A). 



(14) 
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FIG. 3. Comparison of the velocity relaxation rate measured 
in [24] with theoretical predictions using HNC and screened 
Coulomb effective potentials. For a comparison of the exper- 
iment to other theories, see Fig. 5 of |24| . 



FIG. 4. A comparison of classical MD simulations and theo- 
retical predictions for the generalized Coulomb logarithm in 
like-charge electron-ion thermal relaxation. 



where 

For a OCP s = s' . 

Figure [2] shows a comparison between the self-difFusion 
coefficient obtained from using MD particle data in 
the Green-Kubo relations (using the method and code 
described in and predictions obtained using the 

MD, HNC and screened Coulomb effective potentials in 
Eqs. The fi gure shows excellent agreement for 

screened Coulomb in the weakly coupled limit, but this 
breaks down for F > 1, as expected. Likewise, the HNC 
result is excellent in the weakly coupled regime, but this 
agreement also extends into the strongly coupled regime. 
This confirmation that the binary collision picture can 
be extended into the strongly coupled regime through the 
use of an effective potential that accounts for correlations 
is a primary result of this Letter. It provides a prac- 
tical means for evaluating transport coefficients across 
coupling regimes. The results begin to diverge near the 
cross-over to the liquid regime |26| . The curve obtained 
from using the MD extracted effective potential shows 
that only a small part of the disagreement between the 
direct MD data and effective potential theory comes from 
inadequacies of the HNC approximation. 

Comparison with an experiment: Recently Bannasch 
et al |24j measured the velocity relaxation rate in a 
strongly coupled plasma. This experiment was conducted 
in an ultracold neutral plasma held in a magneto-optical 
trap and formed by photoionizing laser-cooled strontium 
atoms. Initially skewed velocity distributions for two spin 
states were formed using optical pumping, and the sub- 
sequent relaxation rate of the two distributions measured 
using laser induced fluorescence j24j . The ions in the sys- 



tem were in a strongly coupled regime, and the electrons 
formed a weakly coupled neutralizing background, pro- 
viding a plasma in which the ion component is, to a good 
approximation, a classical OCP. 

Bannasch et al extracted an average relaxation rate 
(7) from the time-resolved LIF measurements, and com- 
pared the results with predictions of previous theories 
that have the form, 7/0;^ — 0.46r'^/^S, where S is a gen- 
eralized Coulomb logarithm. Figure|3]show a comparison 
of their data with theoretical predictions obtained using 
the HNC and screened Coulomb effective potentials to 
calculate S'^^'^^ The figure shows that, within the mea- 
surement error, the effective potential theory agrees with 
the experiment when correlations are accounted for. 

e"*" — i"*" temperature equilibration: Figure [4] shows a 
comparison of the theoretical predictions and MD results 
for the like-charge electron-ion temperature relaxation 
rate. The generalized Coulomb logarithm is shown in 
the figure, which is inferred from dTf,/dt = 2Q^~^ /SUf, 
where Qe-i = —i'meinf,Vf.i{Tf. — Ti)/mi is the energy ex- 
change density for Maxwellian distributions |3D] and S 
enters through v^i from Eq. ([9|. Details of the MD sim- 
ulations and analysis are explained in |25j . The figure 
shows similar accuracy of the effective potential calcula- 
tion for temperature relaxation across coupling regimes 
as was found for the self diffusion coefficient of an OCP. 

Figures [2]-[4] show that the effective potential theory ac- 
curately predicts both experimental and ah initio simula- 
tion data for a variety of different transport coefficients. 
This demonstrates both the flexibility of this approach 
(because it is compatible with the Chapman- Enskog for- 
malism) , and that the binary collision approximation can 
be extended into the strong coupling regime through the 
use of an effective potential that includes correlation ef- 
fects. The approximation was shown to break down at 
very large F, where there is a known transition to liquid 
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behavior and caging effects turn on [25]. Although the 
range of T values for which this approximation is accu- 
rate is sufficient for many plasma physics applications, 
further refinements to the theory can also be envisioned. 
A potentially significant extension would be to account 
for a dynamic response function in both the closure that 
determines and the relationship between the pair 

correlation function and the effective potential, Eq. 
This may provide an effective potential that accounts for 
relative particle velocities (u), and associated wake ef- 
fects. 
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